Embedding the 125 GeV Higgs boson measured at the LHC in an effective MSSM: 
possible imphcations for neutrahno dark matter 
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We analyze the phenomenological consequences of assuming that the 125 GeV boson measured at 
the LHC coincides with one of the two CP-even Higgs bosons of an effective Minimal Supersymmetric 
extension of the Standard Model at the electroweak scale. We consider the two ensuing scenarios 
and discuss critically the role of the various experimental data (mainly obtained at colliders and 
at B-factories) which provide actual or potential constraints to supersymmetric properties. Within 
these scenarios, properties of neutralinos as dark matter particles are analyzed from the point of 
view of their cosmological abundance and rates for direct and indirect detections. 



I. INTRODUCTION 

It is remarkable that the neutral boson with a mass 
of 125-126 GeV, measured at the LHC in the dipho- 
ton, ZZ, WW and tt channels (hereafter denoted by 
Hi2b) [l|, can be interpreted as the Higgs particle of the 
Standard Model (SM) . However, due to the well known 
problems of quadratic divergences related to the Higgs 
mass, a pressing question is whether this newly discov- 
ered Higgs-like particle can be interpreted within a su- 
persymmetric extension of the SM, where the problem 
of divergences would be solved by boson-fermion loop 
cancellations. Should this be the case, a very rich and 
intriguing phenomenology would open up J2l-lld| . 

Here we investigate this possibility in detail, also in 
connection with possible implications for supersymmetric 
candidates of dark matter (DM) in the Universe. We em- 
ploy a simple supersymmetric model, which we already 
used in previous analyses |8l. ll7Hl9{ . consisting in an effec- 
tive Minimal Supersymmetric extension of the Standard 
Model (MSSM) at the electroweak (EW) scale, where the 
usual hypothesis of gaugino-mass unification at the scale 
of Grand Unification (GUT) of a supcrgravity (SUGRA) 
model, is removed; this effective MSSM is very manage- 
able, since expressible in terms of a limited number of 
independent parameters. 

The Higgs sector of this MSSM has two Higgs doublets, 
which generate, by spontaneous symmetry breaking, two 
vev's: vi and W2- These provide masses to the down-type 
quarks and the up-type quarks, respectively. As usual, 
an angle /3 is introduced and defined as tan/3 = v^jvx. 
This Higgs sector contains three neutral bosons: two CP- 
even, h, H, and a CP-odd one A. The two CP-even Higgs 
bosons are defined, in terms of the neutral components of 



° ^ sin aH^, 



the original Higgs doublets, as: H = cos aH^ 
h = — sin aHi + cos aif2 • In the diagonalization of their 
mass matrix, the mass hierarchy rrih < mn is imposed, 
and the angle a is taken in the range [— 7r/2,7r/2]. 

The lower bound on mh can be established using the 
LEP data on the search for Higgs particles [20[. Con- 
trary to the usual assumption (employed in most of the 
literature until very recently) of assuming for this lower 
limit the Standard Model bound rrih > 114 GeV, in our 
previous works (and in the present one) we take the ac- 
tual LEP constraint on the Higgs-production cross sec- 
tions, that can be translated into a bound on the quan- 
tity sin (a — /3) as a function of nih (this quantity repre- 
sents the ratio of the cross section for the Higgs-strahlung 
process e+e" -^ Zh to the corresponding SM cross sec- 
tion; a complementary bound arises from the process 
ZA, and it is also taken into account). In Refs. 
9| we showed that the LEP limit actually allows 
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the mass tua to reach values as low as 90 GeV in some re- 
gions of the supersymmetric parameter space in MSSM, 
where the phenomenology of ncutralino DM is particu- 
larly interesting. 

The ATLAS and CMS data exclude that the boson A 
can be identified with the new particle H125, but we are 
left with the two options: either H = H125 (hereafter 
denoted as scenario I) or h = i/125 (scenario II). 

As pointed out in Ref . [8[ , scenario I arises naturally in 
the supersymmetric scheme considered in Refs. |17Hl9l |. 
when m/i is taken as light as possible (compatibly with 
the mentioned LEP bound). In fact, in this regime one 
has mhr~.rnA^ (90 - 100) GeV, and tuh ^ (115 - 130) 
GeV [81 . This scenario has also been discussed in Refs. 

m\ 



13l . |15| . As remarked in Ref. 



12, 21 



, the light 

h boson of this scenario could be the origin of the small 
excess of Higgs-like events observed at LEP [22 1. 



The second option is represented by scenario II: this 
scenario occurs, when the Higgs-like boson observed 
at LHC is identified with the lighter CP-even boson 
h within the MSSM [1-0, SO, U^^- In this case, 
rrih — (125 — 126) GeV and A, H can also decouple sub- 
stantially from /i, but with rriA — rnu- 

In this paper we analyse separately scenarios I and II, 
critically discussing the role of the various experimental 
data (mainly obtained at colliders and at B-factories) 
which provide actual or potential constraints to super- 
symmetric properties. For each scenario the properties 
of neutralinos as DM particles are then analyzed from 
the point of view of their cosmological abundance and 
rates for direct and indirect detections. 

The scheme of the presentation is the following. In 
Sect. HI] the features of the employed MSSM are de- 
scribed, in Sect. IIIII a full list of conceivable constraints 
on the model is introduced. Results are given in Sects. 
llVIVIVII and conclusions in Sect. EIIl 

II. EFFECTIVE MSSM 

The supersymmetric model we consider here is an effec- 
tive MSSM scheme at the electroweak scale, with the fol- 
lowing independent parameters: Mi, M2, Afs, /x, tan/3, 
ruA, m^,2, mj, TO[^^^ , "^Fi^^, ™fi., "^fn and A. We 
stress that the parameters are defined at the EW scale. 
Notations are as follows: Mi, M2 and M3 are the U(l), 
SU(2) and SU(3) gaugino masses (these parameters are 
taken here to be positive), ^ is the Higgs mixing mass 
parameter, tan/3 the ratio of the two Higgs v.e.v.'s, wla 
the mass of the CP-odd neutral Higgs boson, rriq^^ is a 
squark soft-mass common to the squarks of the first two 
families, mi is the squark soft-mass for the third family, 
mj and mj are the slepton soft-masses common 
to the L,R components of the sleptons of the first two 
families, m^^ and TTif^ are the slepton soft -masses of the 
L,R components of the slepton of the third family, A is a 
common dimensionless trilinear parameter for the third 



family, Aj^ = A^ = Ami and A^ = A(r 



mf„)/2 (the 



trilinear parameters for the other families being set equal 
to zero). In our model, no gaugino mass unification at a 
Grand Unified scale is assumed, and therefore Mi can be 
sizeably lighter than M2. Notice that the present version 
of this framework represents an extension of the model 



discussed in our previous papers |17H19| . where a com- 
mon squark and the slepton soft mass was employed for 
the 3 families. 

The linear superposition of bino S, wino W^^' and of 
the two Higgsino states H^, H2 which defines the neu- 
tralino state of lowest mass m^ is written here as: 



We assume R-parity conservation to guarantee that the 
lightest supersymmetric particle is stable (we consider 
only models where this is the neutralino). 

Within our model we calculate all the quantities nec- 
essary to impose the constraints discussed in Sect. IIIII 
and the cross sections relevant for direct and indirect de- 
tection of DM neutralinos: the neutralino-nucleon cross 
and the thermally averaged product of 



(nucleon) 
scalar 



section a^ 

the neutralino pair annihilation cross section times the 

relative velocity (o'ann'y)- 

The neutralino-nucleon scattering takes contributions 
from (/i, A, H) Higgs boson exchange in the t-channel 
and from the squark exchange in the s-channel; the A- 
exchange contribution is suppressed by kinematic effects. 
This cross section is evaluated here according to the for- 
mulae given in Ref. [23| . For the crucial coupling pa- 



rameter gd entering the Higgs boson exchange ampli- 
tude, we take its reference value qd.tei = 290 MeV em- 
ployed in our previous papers 



18, a 



We recall that 
this quantity is affected by large uncertainties 1241 with 




(5d,max/ffd,rGf)^ = 3.0 and {gd,min/gd,rcff = 0.12 PJ, Hi 

We also calculate (crannw)int which is the integral of 
(cannw) from the present temperature up to the freeze- 
out temperature T^, since this quantity enters the neu- 
tralino relic abundance (and, for dominant s-wave anni- 
hilation, implies (crannW)int = 2^/ (CannW)) : 



n^h^ 



^f 



9.9 • 10-28 cm^s-i 



gi,{xf) (fannW)int 



(2) 



where Xf is defined as Xf 



^/Tf and gi,{xf) denotes 



a3H°+aiH°. 



(1) 



the relativistic degrees of freedom of the thermodynamic 
bath at Xf. 

The values of VL^h^ . as obtained from Eq. ([2]), are em- 
ployed to exclude neutralino configurations which would 
provide values exceeding the upper bound for cold dark 
matter (CDM), {^cDMh'^)inayii and to rescale the local 
neutralino density p-^, when Q,-^h? turns out to be below 
the lower bound for CDM, {^CDMh'^)inin- In the latter 
case, we rescale p^ by the factor ^ = Pxl POi where po is 
the total local DM density ; S, is conveniently taken as 
e = min{l,r!x/iV(^CDA//i2)mi„} 0. 

In the present analysis, (^cdm ^^)min and 
(ricnA/^2)max are assigned the values: {^cDMh?)ram = 
0.11, (ric'_Dj\//j^)max = 0.13 to conform to the new 
measurements by the Planck Collaboration [26| . 

Since the rates of DM direct detection and those of 
processes due to pair annihilation (with the exclusion of 
processes taking place in macroscopic bodies) are propor- 
tional to ppf and p^ , respectively, in the following we will 

consider the quantities ^'''scalar'™ and ^^(cranni')- 

We calculate Higgs-boson masses and production cross 
sections using FcynHiggs [27 1. 



III. CONSTRAINTS 

We give here a listing of requirements and constraints 
derived from a rich set of experimental data. In Sects. 
IIIIAI - llII Bl are indicated the requirements which are es- 
sential to qualify the model we are considering. Sect. 
unci reports other constraints which can potentially 
bound the physical region of the supersymmetric pa- 
rameter space, but whose implications for our model are 
more involved and thus possibly less compelling; we will 
also explicitly consider the possibility of relaxing some of 
them. This aspect will be discussed later on. 



A. Constraints from the CERN e+e" collider LEP2 

These constraints take into account all data on super- 
symmetric and Higgs boson searches |28[ done at LEP2 
(some of which are improved by those obtained at the 
Tevatron and LHC, as discussed in the next subsection) 
as well as the upper bound on the invisible width for the 
decay of the Z-boson into non Standard Model particles: 

r(z -^ xx) < 3 Mev iili^. 



Bounds from search for Higgs decaying to tau 
pairs. An upper bound in the plane tua — tan/3 is ob- 
tained, in our model, in an indirect way from the data 
reported by the CMS Collaboration in Ref. [3^. A con- 
sistency check of our procedure has been performed, by 
using the upper bound on the production cross section 
reported in Ref. 33[ to obtain the corresponding upper 
bound in the plane m^ — tan/3. 

Bounds on squark masses of the first two fam- 
ines and on the sbottom mass. These bounds are 
taken from the CMS official analysis of Ref. 34 1 . 

Bounds on the stop mass. These bounds are taken 
from the official ATLAS analyses of Ref. 35| for heavy 
stops and of Ref. [36| for light stops. 

Decay Bg — > /i^ + /i^ . We implement the constraint 
recently derived by the LHCb Collaboration in Ref. [37| : 

1.1 X 10"'' < BR{Bs -^ f-i+ + fi-) < 6.4 X 10"^ (8) 

This is a 95 % C.L. hmit. 

Search for the decay t ^> b + H+ . 

Whenever relevant, we have adopted the ATLAS 2-cr 
upper bound on the branching ratio B(i — )• b + H^) as 
reported in Ref. 38 1. 



B. Constraints from the Tevatron and the LHC 



Bounds on searches for supersymmetry from Tevatron 
and LHC are implemented as schematically outlined be- 
low. The observation of a Higgs-like particle seen at 
the LHC imposes specific requirements on the signal 
strengths factors for the production and decay of this 
boson, which have been applied as discussed below. 

Signal strength factors for Higgs produc- 
tion/decay, hi the spirit of the present analysis, in 
the scanning of the supersymmetric parameter space we 
select configurations which satisfy the following require- 
ments, as established by the most recent results at LHC 
El: 



0.61 < R^^ < 1.57 

0.75 < Rzz < 1.47 

0.44 < Rww < 1.24 

0.21 < Rrr < 1.90, 



where the ratio R^j is defined as: 



R-^ 



a{p + p^ Hi25)BR{Hi25 ^7 + 7) 

crsAliP + P -> Hi2;^)BRsm{Hi25 -5-7 + 7)' 



(3) 
(4) 
(5) 
(6) 



(7) 



and similarly for the other final states. Notice that the 
ranges of Eqs. ((SHS]) are 2a intervals. 



C. Constraints from B factories and from (g - 2)^ 
measurements 

Flavor physics experiments are providing stringent 
bounds on many physical processes that can be sizably 
affected by supersymmetric virtual corrections. Here we 
list the most relevant ones for our analysis in the specific 
model we are assuming. 

Measurement of the branching ratio oi h ^ s+j. 

The rate for the branching ratio of the process 5 — > ,s + 7 
is taken here as 2.89 x lO"-* < BR{b -^ 57) < 4.21 x 10"''. 
This interval is larger by 25% with respect to the exper- 
imental determination 39[ in order to take into account 
theoretical uncertainties in the supersymmetric contri- 
butions [40[ to the branching ratio of the process. For 



the Standard Modd calculation, we employ the NNLO 
results from Ref. 



41| 



Search for the decay S — > r -f j/. We use here the 
range 0.38 x 10"'' < BR(B ^ t + i') < 1.42 x lO"'' (world 
average at 95 % C.L.) [ii ]. 

Search for the decay B ^ D + t + ly. A new range 
for the quantity R{D) = BR{B -^ Dtv)/BR[B -^ Dev) 
has been established by the BABAR Collaboration [43J : 
30.0 X 10-2 < R{D) < 58.8 x lO'^ (2cr interval). 

Muon anomalous magnetic moment (5 2)^. We 

take the conservative 2ct range 3.1 x 10"^° < Aa^ < 



Scenario I 


tan/3 


(4,6) 


M 


(1800, 2000) GeV 


Ah 


(40, 80) GeV 


Ah 


(180, 800) GeV 


Ah 


~ 2000 GeV 


"l912 


(1400, 1600) GeV 


mj 


(1400, 1600) GeV 


m,- , m,- 

'12,1,' '12, R 


~ 500 GeV 


mf ^ , nifj, 


(120, 200) GeV 


rriA 


(100, 120) GeV 


\A\ 


(2.5, 2.8) 



TABLE I: Values and intervals for the MSSM parameters out- 
lined by the LHC bounds for Scenario I. 



"^'^ of the 



47.9 X 10~i° for the deviation Aa^ = 
experimental world avera ge o f a^ = (g^ — 2)/2 4^ from 
the theoretical evaluation [45J (in the latter we estimate 
the leading hadronic vacuum polarization contribution 
in the Standard Model by combining the two determi- 
nations estimated from e^e^ and r-dccay data). We 
evaluate the supersymmctric contributions to the muon 
anomalous magnetic moment within the MSSM by using 
the formulae in Ref. 



IV. SELECTION OF SUPERSYMMETRIC 

CONFIGURATIONS AND NEUTRALINO RELIC 

ABUNDANCE 



Here we provide the results of our analysis for scenar- 
ios I and II, regarding the selection of supersymmctric 
configurations and the ncutralino relic abundance. The 
mass interval for the LHC Higgs-like particle is taken 
here as 123 GcV < m^i^s < 129 GeV. 



A. Scenario I: J/ = Hri5 

This scenario is defined by identifying the heavier CP- 
cvcn Higgs neutral boson H of the MSSM with the LHC 
Higgs-like particle, i.e. H = H125. This implies that 
the mass interval for H125 obtained at the LHC has to 
be attained by H (123 GcV < m^ < 129 GeV), and 
this entails to/j ^ itia — (90 — 100) GcV, as already 
noticed above. The LHC constraints on the production 
rates in the various channels detailed in Eqs. 1^^ select 
the sector of supersymmctric parameter space reported 
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FIG. 1: Scenario I - Signal strength factors as defined in Eq. 
(O for the production and decay of the heavy Higgs H when 
123 GeV< niH < 129 GeV). (a): R-y--,; (b): Rzz\ (c): Rww\ 
(d): Rtt- The horizontal lines denote the allowed intervals 
obtained at the LHC, and given in Ea. (|3l4l5l6[) . 



in Table n ^ 

The most peculiar feature of this region is represented 
by the high values of the parameter /x, a property which 
agrees with the findings of Refs. 



Il3lll5| and appears to 
be related to the constraint imposed by Eq. ^ on R-y^p 



as remarked in Ref. 131. We also notice that the sector 



of parameter space defined in Table U has some similar- 
ities with the scenario denoted by low-Mjj in Ref. [15J . 
though it differs in one important feature: in our case 
the slepton masses (most notably the mass parameters 
for tl and f/j) are significantly lower. These are the pre- 
requisites for having configurations where the neutralino 
relic abundance does not exceed the cosmological bound. 

The properties of the solutions we have found are dis- 
played in Figs. [T] - m In Fig. [T] we show the signal 
strength factors for Higgs production and decay at the 
LHC; we notice how a sizable subset of our population 
of supersymmetric configurations fits quite well all LHC 
data on these factors. This population satisfies also the 



For the lower bound on the slepton masses we use here the LEP 
values mj- > 80-100 GeV (depending on flavour) [3011 . These 
lower bounds actually depend on the condition that m^—my^^ > 
0(3— 15) GeV. If these conditions are not met, the slepton lower 
bound can decrease to about 40 GeV, with relevant implications 
for the neutralino phenomenology, as discussed in Ref. [4711 . 
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FIG. 2: Scenario I - Some of the experimental constraints dis- 
cussed in Section Uni are compared to the corresponding the- 
oretical expectations for the supersymmetric configurations 
reported in Table |T] and Fig. [T] Panels (a)-(c) correspond 
to collider constraints: (a) LEP bound on the higgs produc- 
tion cross-section, reported in terms of the coupling factor 
sin^(Q — p) [2g|; (b) CMS bound on Higgs production and 
subsequent decay into rf [33| (the production cross section 
refers Xq (^ = A unless tua — rn^ or niA — rnn, in which 
case (f) — A,h OT <f) = A,H, respectively); (c) ATLAS upper 
bound on the branching ratio BR(i — > bH^) [38[. Panel (d) 
shows the extent of deviations from the two constraints on 
BR(fe^s + 7) and (g-2)^. 



other relevant constraints from colliders as depicted in 
panels (a)-(c) of FigH] Panel (d) of Fig. [51 instead, 
shows that predictions for BR(6 — > s -I- 7) and {g — 2)^ 
in scenario I deviate from the experimental bounds dis- 
cussed in Sect. IIII CI A minimal deviation occurs for 
(g — 2)^, whereas a deviation of about 4cr occurs for 
BR(6 -^ s -I- 7). Therefore scenario I, which is perfectly 
viable as far as accelerators data are concerned, is in 
tension with the experimental bounds when also flavor 
physics determinations are included (this will not be the 
case for scenario II, as discussed below). Contrary to ac- 
celerator physics constraints, these are indirect bounds 
and rely to some degree of cancellation of various terms 
[19|, which may not be fully under theoretical control. 
We therefore discuss the implications of scenario I for 
dark matter, nevertheless reminding that this scenario 
exhibit a significant level of tension with indirect bounds 
on supersymmetry. 

In the present scenario the neutralino mass sits in 
the range: m^ ~ (40 — 85) GeV. As shown in Fig. |31 





Scenario 


II 




tan^ 


(4, 20) 




M 


(100, 400) GeV 




Ml 


(40, 170) GeV 




M2 


(100, 1000) GeV 




Ms 


~ 2000 GeV 




rug,^ 


(700, 2000) GeV 




nil 


(700, 1200) GeV 


"^fl2,^ 


"^h2.R^'^fL,rnf^ 


(80, 1000) GeV 




rriA 


(200, 1000) GeV 




1^1 


(1.5, 2.5) 



TABLE II: Values and intervals for the MSSM parameters 
outlined by the LHC bounds for Scenario II. 



most of our configurations have a sizable neutralino relic 
abundance. Fig. |3] iUustrates the contributions of differ- 
ent annihilation channels to the integrated cross section 
(o'annw)mt; wc noticc that, as anticipated, light sleptons 
are instrumental in keeping the annihilation cross sec- 
tion large enough to comply with the experimental up- 
per bound on O-^/i^, since diagrams with exchange of a 
slepton dominate (iTannf) mt over the whole range of the 
available neutralino masses, with the exception of a small 
mass range around my^ ~ mA/'2-, where resonant annihila- 
tion through A exchange can become important. On the 
other hand, Z~boson exchange remains sub-dominant 
even close to the pole in the corresponding annihilation 
cross section, m^ ~ Mz/'^, since the Z boson couples 
to the neutralino only through its Higgsino components, 
while in this scenario the neutralino is a Bino of ex- 
tremely high purity, due to the very large values required 
for the /x parameter, as specified in Table HI Finally, for 
the same set of supersymmetric configurations we show 
in Fig. [5]the ratios [(crannw)i/(a'annf)tot]T=o between the 
neutralino annihilation cross sections times velocity to 
the final states i = tt, bb and the total annihilation cross 
section times velocity, both calculated at zero tempera- 
ture. The latter quantities are relevant for the evaluation 
of indirect signals, as discussed in Section IVll As shown 
in the plot, annihilation to rf (driven by the exchange 
of light staus) is dominant, with a sub-dominant contri- 
bution from the bb annihilation channel (which, as in the 
case of ((Ja_nnv)inti Can become sizeable through resonant 
annihilation through A exchange). 



B. Scenario II: h = H125 

This scenario is defined by the alternative choice, that 
identifies the lighter CP-even Higgs neutral boson h 
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FIG. 3: Scenario I - Neutralino relic abundance as a func- 
tion of the neutralino mass for the supersymmetric configu- 
rations reported in Table |T] and Fig[l] The horizontal line 
represents the upper bound (f2ci3M^^)max = 0.13 from the 
Planck Collaboration [2g| on the cold dark matter content in 
the Universe. 



FIG. 4: Scenario I - Fractional contributions of different an- 
nihilation channels to the integrated neutralino cross section 
times velocity (aannv) int for the supersymmetric configura- 
tions reported in Table [T] and Fig[l] (Green) dots: XX ~^ ff 
through slepton exchange; (red) crosses: xX ~^ ff through 
Higgs exchange; (blue) open circles: xX ^^ ff through Z- 
boson exchange. 



of the MSSM with the LHC Higgs-like particle, i.e. 
h = H125. This therefore implies: 123 GeV < nth < 
129 GcV. The scan of the MSSM parameter space that 
produces a population of configurations satisfying all re- 
quirements and constraints mentioned in Sect. |llT] iden- 
tifies the sector outlined in Table |IT1 The features of this 
population are displayed in Figs. [51- 1101 

Fig. [5] shows how the requirements for the signal 
strength factors are verified for our configurations. The 
constraint derived from LHC searches for Higgs decaying 
to a tau pair implies for the mass of the CP-odd Higgs A 
the lower bound: to^ > 300 GeV, as indicated by panel 
(a) of Fig. [71 Fig. [Tljb) shows that, at variance with the 
previous case, in scenario 11 the constraints BR(5 — > s+7) 
and {g — 2)^ are satisfied. It is worth noting that also 
the bound on the branching ratio for the invisible de- 
cay h 



X + X 



not explicitly discussed before, is 
respected. 

The plot of Fig. [51 displaying the neutralino relic abun- 
dance versus the neutralino mass, shows that m^ has the 
lower limit m^ > 30 GeV and that there exists a break in 
the range 70 GeV < rriy^ < 85 GeV, this interval being dis- 
allowed by the requirement that fl^Ii^ < {i^cDMh^)ma.x- 
In turn, this property is due to the strong enhancement 
in the pair annihilation amplitude when m^ runs over 
the values ^rrifi, \mz. 



Fig. [21 displays the various contributions to (crann*^)mt 
and shows that dominances in the annihilation ampli- 
tude are as follows: a) dominance of annihilation to 
fermions through Z-exchange in the range 30 GeV < 
TO^ < 60 GcV; b) dominance of annihilation to fermions 



through light scalar Higgs-exchange for m^ 



mhl2 



; c) dominance of annihilation to W'^W" for m^ > 
mw . To conclude the discussion, Fig[TOl shows the ra- 
tios [(cranni')i/(o'anni')tot]T=o between the neutralino an- 
nihilation cross section times velocity to the final states 
i = TT,bb,W~^W' ,ZZ,Zh and the total annihilation 
cross section times velocity, both calculated at zero tem- 
perature: as shown in the plot, rf dominates when m^ < 
65 GeV, W'^W~ prevails when myy < m^ < mz + rrih^ 
and Zh dominates at larger masses. 



DIRECT DETECTION. 



We turn now to the evaluation of the relevant quantity 
for DM direct detection, C^'scaia™" • "^^"^ values for this 
quantity are shown in the scatter plot of Fig. II II together 
with the regions pertaining to the signals measured by 
the experiments of DM direct detection of Refs. [5]J, I54I 
(other experimental results showing an excess of events 
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m^(GeV) 

FIG. 5: Scenario I - Ratios [(crann«)i/(''"ann'')tot]T=o between 
the neutralino annihilation cross sections times velocity to the 
final states i = tt, bb and the total annihilation cross section 
times velocity, both calculated at zero temperature, for the 
supersymmetric configurations reported in Tableland Fig[T] 
(Red) dots: tt final state; (blue) crosses: bb final state. 



compatible with a positive signal are reported in Refs. 



53|,l54[). In particular, in Fig. [TT](red) crosses represent 



configurations found in the set of Scenario I, while (blue) 
dots correspond to configurations found in Scenario II. 
The experimental domains shown here were obtained by 
using for the velocity distribution function of the galac- 
tic dark matter those pertaining to a standard isothermal 
sphere with po = 0.30 GeV cm~'^, vq = 220 km sec""'^, 
with Uesc = 650 km sec"! for the DAMA/LIBRA ex- 
periment and Ucsc = 544 km scc^^ for CRESST and for 
specific sets of experimental parameters (quenching fac- 



5lll5^- Including 



tors and others), as discussed in Refs. 

uncertainties of various origin, the experimental regions 

would expand as indicated for instance in Fig. 7 of Ref. 

0. 

One notices that the set of configurations found in the 
scenario I generate very low rates for direct detection of 
relic neutralinos. Thus in this scheme neutralinos does 
not appear be responsible for the signals meas ured by the 
experiments of DM direct detection of Refs. 



51-54 



It is worth stressing that these conclusions rest heavily 
on the results recently obtained from colliders; in par- 
ticular, very constraining are the conditions expressed in 
Eqs. PI4I5I6P and the bounds implied by the search for 
Higgs decay into tan pairs, that constrain the parameter 
/i to be very large and tan/3 small. Should these con- 
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FIG. 6: Scenario II - The same as in Fig. [T] except that here 
123 GeV< ruh < 129 GeV. 



straints significantly relax in the future, as a consequence 
of further experimental data and analyses from colliders, 
the theoretical values of ^cr^"^;^^"" would compare to the 
data of DM direct detection much more favorably, as for 
instance depicted in Fig. 5 of Ref. [8|. 

In the case of scenario II, in view of the experimen- 
tal uncertainties mentioned above and of the theoretical 
uncertainties related to the parameter gd (see Sect. [Ill j 
the gap between the experimental regions and the scat- 
ter plot shown in Fig. [TT] could somewhat narrow down. 
Most of the theoretical values shown in Fig. [TT] are in 
tension with the experimental bounds given by other DM 
experiments (see for instance Ref. 55ll56|). 



VI. INDIRECT DETECTION. 

In order to study the capability of indirect signals to 
probe neutralino dark matter in scenario I and scenario 
II, we discuss the exotic component in cosmic rays rep- 
resented by antiprotons, and the contribution to the so- 
called isotropic gamma-ray background (IGRB) due to 
the production of gamma-rays at high latitudes from an- 
nihilation in our Galaxy. 

Antiprotons are potentially able to provide quite 
strong bounds on dark matter annihilation in our Galaxy. 
We therefore calculate the antiproton production in both 
scenarios and compare them with the PAMELA mea- 
surements of the absolute antiproton flux [57j . Similar 
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FIG. 7: Scenario II - Two of the experimental constraints 
discussed in Section IIIII are compared to the corresponding 
theoretical expectations for the supersymmtric configurations 
reported in Table HTl and in Fig. |6l (a) CMS bound on Higgs 
production and subsequent decay into tt [321 (the production 
cross section refers to (j) = A unless tua — rrih or rriA — rnn, 
in which case (j) = A,h or (f) = A,H, respectively); (b) the 
two constraints on BR(& — ?> s + 7) and (g — 2)^, which in this 
scenario are simultaneously satisfied. 



bounds can be obtained with the BESS-Polar II deter- 
mination 58|. Fig. [T^ shows the antiproton fluxes in the 
first PAMELA energy bin {Tp = 0.28 GeV) for the config- 
urations of scenario I (red crosses) and scenario II (blue 
circles). The left panel refers to a galactic propagation 
model with the MED values of propagation parameters 
[59| : the right panel refers to the MAX set of parame- 
ters [59[. The MAX set refers to the configuration in 
the space of propagation parameters which provides the 
largest antiproton fluxes (mostly due to a large volume of 
the cosmic-rays confinement region) , while being allowed 
by B/C measurements [59[. 

The upper long-dashed line denotes the 95% C.L. 
bound by using the PAMELA data 57[ and adding in 
quadrature a 40% theoretical error on the theoretical de- 
termination of the antiproton background. This generous 
allowance is taken under consideration because of uncer- 
tainties in the nuclear cross sections relevant for the sec- 
ondary production [59| . The modification of the bound 
when a smaller estimate of the theoretical uncertainty 
(20%) is adopted |59[ is shown by the short-dashed line. 
We notice that antiprotons are far from bounding the 
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FIG. 8: Scenario II - The same as in Fig[3] for the super- 
symmetric configurations reported in Table |II] and in Fig. |H] 
(Scenario II in the text). 



configurations of both scenario I and II. This is due to 
the fact that the dominant channel of annihilation in a 
large portion of the parameter space of both scenarios is 
a leptonic one (namely, tt) which is unable to produce 
a relevant amount of antiprotons. Only those configu- 
rations of scenario I where the bb final state dominates 
(very few configurations with a neutralino mass close to 
55 GeV, as seen in Fig. [5]) and the configurations of 
scenario II where the gauge-bosons final state (for neu- 
tralino masses above 80 GeV, as seen in Fig. [TU|) domi- 
nates are able to produce an antiproton flux that reaches 
its maximal values. Dominant hadronic (bb) final states 
for neutralino masses below 70 GeV are accompanied by 
small values of the neutralino relic abundance: this has a 
strong impact in reducing the antiproton flux, due to the 
squared appearance of the rescaling factor ^ in indirect 
signals (since they depend on p^). 

Current antiprotons bounds therefore do not con- 
strain our supersymmetric conflguration, neither for the 
MIN nor for the MAX set of propagation parameters. 
Prospects for future searches are shown in Fig. [T^by the 
dotted line, which refers to an expected reach of AMS 
j60| . We estimated AMS capabilities by taking into con- 
sideration the following facts: AMS data on antiprotons 
will likely reach a level of a few percent uncertainty; AMS 
will determine the fluxes of cosmic rays species to an un- 
precedented level, and this will help in reducing also the 
theoretical modeling of galactic cosmic rays propagation. 
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FIG. 9: Scenario II - The same as in Fig. [J] for the su- 
persymmetric configurations reported in Table [III and in Fig. 
[B] (Green) dots: XX ~^ ff through slepton exchange; (red) 
crosses: xX ~^ ff through Higgs exchange; (blue) open cir- 
cles: XX ~^ ff through Z-boson exchange; (dark green) open 
squares: xX -^ WW; (purple) open triangles: xX ~^ ^^^ 
(cyan) filled circles: xX ~^ ^h; (grey) filled triangles: xX ^^ 
hh. 



Determination of the boron-to-carbon (B/C) ratio will 
be especially relevant. By considering a total (theoreti- 
cal -|- experimental) uncertainty on the antiproton fluxes 
after AMS, we can estimate a bound (in case of non ob- 
servation of deviation from the expected background) at 
the level of the dotted lines in Fig. [12) this would al- 
low to probe a fraction of the parameter space, both for 
scenario I and scenario II, in the case of relatively large 
values of the propagation parameters (right panel of Fig. 
[T2)) . This capability is further illustrated in Fig. [T3l 
where two representative antiproton fluxes (one for sce- 
nario I and one for scenario II) are reported. The two 
fluxes refer to the best-choice occurring in our param- 
eter space, but are nevertheless representative of those 
configurations with fluxes in excess of the AMS reaching 
capabilities shown in the right panel of Fig. [121 Dark 
matter fluxes like those shown in Fig. [T3|will easily rep- 
resent a detectable signal in AMS, considering that they 
produce visible excesses over the background (denoted by 
the solid line, while the dashed lines bracket a ±10% un- 
certainty) in most of the energy spectrum. We also stress 
that AMS will have a very large statistics and therefore 
and excess like those shown in Fig. [T3|will be detected as 
a deviation in a large number of experimental bins, thus 



FIG. 10: Scenario II - The same as in Fig. [5] for the su- 
persymmetric configurations reported in Table HIl and in Fig. 
[B] (Red) dots: rf final state; (blue) crosses: bb final state; 
(dark green) open squares: W'^W~ final state; (purple) open 
triangles: ZZ final state; (cyan) filled circles: Zh final state. 



making the evidence of a signal potentially quite clear. 
The major limitation remains the ability to reduce the 
theoretical uncertainties on the background to a suitable 
level, as discussed above. 

Concerning the indirect signal in terms of gamma-rays, 
Fig. [T3I shows the flux of gamma rays produced by galac- 
tic dark matter annihilation at high latitudes for both 
scenario I and scenario II. The contribution to the IGRB 
has been calculated for an Einasto profile of the dark 
matter density, but different profiles predict only slightly 
different fluxes 



6lj . since we are looking here at high 



galactic latitudes. 

The signal fluxes in both scenario I and scenario II 
are relatively small, when compared to the current up- 
per bounds on the IGRB, obtained by considering the 
Fermi-LAT measurements [631 and the best-fit of var- 



ious contributions to the IGRB [61j: misaligned AGN 
[63j . star-forming galaxies J6J], unresolved milli-second 
pulsars 65 1, BL Lacertae [66| and flat-spectrum radio 
67J |. The upper bound at the 95% C.L. is shown 



quasars 



in Fig. [Ujby the horizontal dashed line. The figure shows 
the Hux at two representative energies, corresponding to 



two different energy bins of the Fermi-LAT analysis [63 : 
the left panel refers to E^ — 1.2 GeV, the right panel to 
E^ = 9.4 GeV. 

We notice that the contribution to the IGRB of as- 
trophysical origin suffers of large uncertainties: in de- 
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FIG. 11: Neutralino-nucleon coherent cross section times 
the rescaling factor 5a"^""j^™"\ (Red) crosses: supersymmet- 
ric configurations plotted in Fig[l] (Scenario I in the text); 
(blue) dots: supersymmetric configurations plotted in FigE] 
(Scenario 11 in the text). The hatched areas denote the 
DAMA/LIBRA annual modulation regions [51(: the (green) 
vertically-hatched region refers to the case where constant 
values of 0.3 and 0.09 are taken for the quenching factors of Na 
and I, respectively [49(]; the (red) crossed-hatched is obtained 
by using the energy-dependent Na and I quenching factors 
as established by the procedure given in Ref. [50|. The gray 
regions are those compatible with the CRESST excess [52l |. 
In all cases a possible channeling effect is not included. The 
halo distribution function used to extract the experimental 
regions is given in the text. For other distribution functions 
see fH 



riving the bounds shown in Fig. [14] we have adopted 
the central-value determinations of the different sources 
of background, as reported in Ref. 



6l| . If (just) some 



of these background fluxes are allowed to fluctuate up 
(especially the recently determined gamma-ray flux orig- 



inating from misaligned AGN 63 1) the ensuing bounds 



can become quite constraining 



61|. 



VII. CONCLUSIONS 

The attempt of interpreting the neutral boson (i?i25) 
measured at the LHC in the diphoton, ZZ, WW and 
TT channels, and with a mass of 125-126 GeV, in terms 
of the effective Minimal Supersymmetric extension of the 
Standard Model defined in Sect. |TT1 has led us to consider 
two possible scenarios: a scenario I, where the boson H125 



is identified with the heavier CP-even neutral boson H 
and scenario II, where the boson H125 is identified with 
the lighter CP-even neutral boson h. 

The supersymmetric parameter space has been anal- 
ysed also in terms of a full set of constraints derived from 
collider experiments, B-factories, and measurements of 
the muon anomalous magnetic moment. The properties 
of the ncutralino as a dark matter constituent has been 
analysed in both scenarios, considering its relic abun- 
dance and direct and indirect detection rates. 

We have found that in scenario I no solution for su- 
persymmetric configurations exists, unless two indirect 
constraints (BR(6 -> s -f 7) and {g — 2)^) are relaxed. 
If these two requirements are not implemented, solu- 
tions with a physical relic abundance are found in a re- 
gion of the supersymmetric parameter space character- 
ized by low values for the stau mass parameters 80 GeV < 
mj , mj , rufj^ , 77if^ < 200 GeV, and high values for 
the /i parameter: fj. > 1.8 TeV. In the region defined in 
Table [I] the neutralino mass turns out to sit in the range 
niy^ ~ (40 — 85) GeV. The set of configurations found 
in the present scenario generate very low rates for direct 
detection of relic neutralinos (the quantity ^o'g""j^^''°" is 
at the level of (cri^^^il^""^ ^ a few xlO'^^ cm^). The 
same occurs for indirect detections signals: only antipro- 
ton searches, under some optimistic assumptions, may be 
able to test scenario I for neutralino masses close to 50 
GeV. For this to be reachable, a somehow large cosmic- 
rays confinement region is required, accompanied by a 
reduction of the total theoretical -I- experimental uncer- 
tainty on the antiproton fiux determination at the level 
of about 10%. AMS [60| is expected to beat this level of 
precision on the antiproton data, and its measurement of 
the fiuxcs of cosmic rays species, especially B/C, could 
help in reducing the uncertainties on the theoretical de- 
termination, allowing to approach the level required to 
study these supersymmetric populations. 

In scenario II we have found a population of config- 
urations which satisfy all requirements and constraints 
mentioned in Sect. IIIIl including the indirect bounds 
coming from BR(6 -^ s + 7) and {g — 2)^. Here the lower 
limit for the ncutralino mass is m-)^ > 30 GeV. The di- 
rect detection rates are shown to be typically rather low; 
though, they could approach the level of the signals mea- 
sured by the experiments of DM direct detection 5ll454 
under special instances for the DM distribution, for ex- 
perimental parameters and/or for significantly large size 
of the neutralino-nucleon coupling. As for the indirect 
signals a situation similar to scenario I occurs: under the 
same, somehow optimistic, assumptions discussed above 
an antiproton signal in AMS may be reachable for ncu- 
tralino masses above 80 GeV. 

A few comments are in order here, regarding the fea- 
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FIG. 12: Antiproton fluxes in the first PAMELA energy bin {Tp — 0.28 GeV) for tfie configurations of scenario I (red crosses) 
and scenario II (blue circles). The upper long-dashed line denotes the 95% C.L. bound by using the PAMELA [57|] data and 
adding in quadrature a 40% theoretical error on the theoretical determination of the antiproton background. The short-dashed 
line shows the same upper bound, with a 20% estimate of the theoretical uncertainty. The dotted line denotes the reaching 
capabilities of AMS [60|, provided the total experimental and theoretical uncertainties are reduced to 10%. The left panel refers 
to a galactic propagation model with the MED values of propagation parameters [59(; the right panel refers to the MAX set of 
parameters [59j . 



tures of the population of relic neutralinos examined in 
the present paper: a) our results apply only to the stan- 
dard situation of thermal decoupling in a standard FRW 
cosmology; in more extended cosmological scenarios, es- 
pecially those with an enhanced expansion rate of the 
Universe the features of these populations are expected 
to be different 681470] ; b) the relic neutralinos considered 
here could constitute only a part of a multicomponent 
DM (another component virould be the one responsible 
for the signals observed until now in DM direct detection 
experiments); c) the derivations presented in the present 
paper rest heavily on the results obtained at colliders: 
many of the analyses pertaining these results are actu- 
ally in progress, thus some of them could be susceptible of 
significant modifications, with the implication of possible 
substantial changes in the our present conclusions. 
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